Variable range hopping in the Coulomb glass 
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We use a mean-field (Hartree-like) approach to study the conductance of a strongly localized 
electron system in two dimensions. We find a crossover between a regime where Coulomb interactions 
modify the conductance significantly to a regime where they are negligible. We show that under 
rather general conditions the conduction obeys a scaling relation which we verify using numerical 
simulations. The use of a Hartree self-consistent approach gives a clear physical picture, and removes 
the ambiguity of the use of single-particle tunneling density-of-states (DOS) in the calculation of the 
conductance. Furthermore, the theory contains interaction-induced correlations between the on site 
energy of the localized states and distances, as well as finite temperature corrections of the DOS. 

PACS numbers: 71.23.Cq, 73.50.-h, 72.20.Ee 



INTRODUCTION 

The unusual temperature dependence of the conduc- 
tance of strongly disordered samples has generated a 
great interest in the past four decades. These systems 
manifest an interplay between disorder and interactions, 
which is ubiquitous in nature, and which renders the 
problem difficult theoretically and rich experimentally 
In these disordered systems, conductance 
occurs via phonon assisted hopping between localized 
states, whose size is roughly the localization length £. 
The measured low-temperature conductance of various 
disordered samples shows temperature dependence dif- 

ferent than simple activation, of the type a ~ e v T > , 
with C < 1- 

Mott's theory ("Mott's law") for Variable Range Hop- 
ping (VRH) @, which ignores the Coulomb interactions, 
gives £ = and To equal to the non-interacting mean 
level spacing in a localization volume 6% = l/Vo£ d ; with 
d = 1,2 or 3 the dimension of the system and uq the 
density-of-states. 

As first pointed out by Pollak Q, interactions affect 
the single-particle tunneling density-of-states (DOS), and 
thus may possibly influence the conductance. Efros and 
Shklovskii (ES) p} obtained an analytic bound for the so- 
called Coulomb gap in the DOS, which was later observed 
experimentally [H, Q. Using the interaction-modified 
DOS to find the conductance (treating the problem as 
essentially non-interacting from this point) one obtains 
("ES's law") C = \i m any dimension d. In ES's law 
To is given by the Coulomb energy of the localized state 
= e 2 /£. The Mott and the ES laws, as well as the 
crossover between the two regimes, were observed in var- 
ious experiments fioj ]. Many qualitative features that 
this complicated system exhibits can be obtained using 
rather simple and heuristic arguments. However, using 
the single-particle DOS in the calculation of the conduc- 
tance, neglecting finite temperature effects on the DOS 
and ignoring correlations between energies and spatial 
locations is not discussed or can not be justified in the 



framework of the heuristic arguments (see 11[ for an 
elaborate discussion). 

In this work we study the crossover between Mott's and 
ES's laws, using the theoretical framework of a Hartree 
theory 12, [l3j]. The Hartree theory takes into account 
interactions in a self-consistent manner. Hence, it gives 
a well defined procedure for calculating various physi- 
cal properties and allows us to develop a clear intuitive 
physical picture. Within the Hartree approximation the 
use of a rcnormalized single-particle DOS is built in, 
with the electrons moving in an effective potential due 
to the other electrons. While various other theories ne- 
glect the interaction-induced correlations between on-site 
energies and position 14, lU , these exist in the theory we 
present: The mean-field equations (l3| contain both the 
electronic site positions and the renormalized energies, 
and therefore correlations between the two exist in the 
self-consistent solution. We emphasize that within the 
Hartree theory the average occupation and renormalized 
energy arc different at each electronic site. Albeit, the 
Hartree method is uncontrolled as there is no small pa- 
rameter in the theory. It is therefore difficult to estimate 
rigorously its limit of validity. Nevertheless, due to the 
long-range nature of the Coulomb interactions, the po- 
tential energy at a given site is renormalized by a large 
number of electronic sites. This suggests that the ap- 
proximation should be valid, and we believe it should be 
applicable for a wide range of system parameters. In- 
deed, the mean-field picture yields the well-established 
results for the Coulomb gap in two dimensions Il3j , and 
captures well the results of aging experiments [161 ] . 

In this work we demonstrate, in two dimensions, that 
the Hartree theory yields the ES's and Mott's laws in- 
cluding the crossover between them, for which we derive 
an analytical scaling relation that holds under generic 
conditions, extending the results of Refs. 17, 18, 1J, 2(ij ]. 
This is a proof that interaction-induced correlations be- 
tween location and energy as well as finite temperature 
corrections are not essential for determining the conduc- 
tance. 

The Hartree approximation involves a numerical pro- 
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ccdurc. Nevertheless, it clarifies our physical understand- 
ing of the complex Coulomb glass. Its success in repro- 
ducing equilibrium (the gap in the DOS) and near equi- 
librium properties (the conductance) and its simplicity 
indicate that it may be usable in the estimates of dy- 
namical properties. Indeed, this was done in [l6j |. 

The structure of the paper is as follows. We first 
present the model, and the heuristic arguments leading to 
the Mott and ES laws. We show the Hartrcc theory gives 
a Coulomb gap consistent with other theoretical and nu- 
merical approaches, and study the temperature depen- 
dence of the DOS. We then explain the mean-field ap- 
proach to determining the conductivity, and the numeri- 
cal procedure involved. Finally, we analyze the crossover 
between the two regimes analytically, under certain as- 
sumptions, and compare this to the numerical results. 
The discussion that follows is centered around the two- 
dimensional case, except for the scaling analysis which is 
done for completeness in arbitrary dimension. 



brackets are replaced by n(\Eij\) . We emphasize that 
the energy difference contains the Coulomb interaction 
with all other charges. The e 2 /r interaction term present 
in the Efros-Shklovskii argument should not be included 
here, as in the master equation approach the charge is 
continuously transferred. Indeed, only without this term 
will the approach lead to detailed balance at equilibrium 
and to the correct result for the Coulomb gap [131 ] , includ- 
ing the coefficient as demonstrated in Fig. [1] Recent 
numerical simulations have confirmed the expected form 
of the three dimensional Coulomb gap [2lj], and it would 
be interesting to see this also for the two dimensional 
case. 



ANALYSIS USING HEURISTIC ARGUMENTS 

The heuristic arguments are based upon the low tem- 
perature approximation that the conductance follows [l[ : 



THE MODEL OF HOPPING BETWEEN 
LOCALIZED STATES 

The model analyzed consists of TV localized states and 
M < N interacting electrons, with a coupling between 
the electrons and a phonon reservoir [l3[ . The states are 
localized at sites r^, i = 1..N, and by assumption the av- 
erage distance between the states is larger than the local- 
ization length length £. The states have different on-site 
energies, ei, due to the disorder in the hosting lattice. For 
simplicity, we will neglect in our analysis fluctuations in 
£ from site to site. Since the states are localized, the elec- 
trons will interact via an unscreened Coulomb potential 
e 2 /\rj — ri\ (to take into account the dielectric constant 
k, one should replace e 2 by — throughout). The disorder 
bandwidth W , will be assumed to be much larger than 
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the Coulomb energy term ^— for the transport calcula- 
tions , with r nn the average nearest-neighbor distance. 
In the self-consistent Hartree approximation the energies 
of the states are renormalized from e, to an energy Ei 
due to the interaction with the other electrons. 

Let us denote the energy difference of the electronic 
system before and after a phonon-induced tunneling 
event by Eij = Ej — Ei, and the distance between their 



-(<-ij/C-|Eijl/T) 



(2) 



locations by 



For weak electron-phonon 



coupling the transition rate 7?„- of an electron from site i 
to site j, is approximately (for Eij < 0)[1|]: 



7° 



\M q \ 2 u p f -f(E J ))e--[l + n(\E. lJ \)i (1) 



where f(E) and n(E) are the Fermi-Dirac and the Bosc- 
Einstein distributions, respectively; M q the electron- 
phonon coupling constant and v v the phononic density- 
of-states. For upward transitions {Eij > 0) the square 



Where 7°- are two optimally coupled sites (i.e, with the 
highest rate between them). To maximize exponential 
functions it is sufficient to minimize the exponent even if 
the prefactors are not optimized. Indeed, it pays some- 
times for the electron to hop a larger distance, thereby 




FIG. 1: A comparison between the Coulomb gap obtained 
by the mean-field approximation, in two dimensions, and the 
theoretical prediction [22] of a linear gap with slope . The 
number of sites used was N=10000, and the Fermi energy was 
Ef = (half filling). The sites were uniformly distributed in 
2 - 

a square, with e 712 = 1, where n is the density. The tem- 
perature was taken to be zero. The on-site disorder was dis- 
tributed uniformly in the interval [— j, 5]. The y axis denotes 
the probability density of the energies Ei. The points are an 
average over 300 instances. 
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finding a state closer in energy and reducing the neces- 
sary inelastic energy transfer [5[ . 

In the Hartree self consistent treatment it is legitimate 
to use the non-interacting concept of DOS and to ask how 
interactions modify it. To determine the shape of the in- 
teracting DOS, v{E), (not to be confused with the non- 
interacting DOS i/q) we assume that the nonintcracting 
and interacting energy-distance relation are equal, i.e., 
for two dimensions E = e 2 /r = \/v(E)r 2 . This yields the 
ES result v(e) ~ E/e 4 . The tunneling DOS is suppressed 
at low energies due to interactions [iSHHEl. No- 
tice that v(E) ~ Vq at Ec = e 4 ^o, an energy above which 
the DOS becomes v$. Hence Ec determines the scale of 
the Coulomb gap. 

Regarding the conductivity of the system, a-priori 
there could be four energy scales involved: the band- 
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width W; the Coulomb energy , with r nn the average 
nearest-neighbor distance; the Coulomb energy over the 

2 -I 

localization length distance JE^ = and 8g = ^j, the 
level spacing in a box of size £. However, it is clear that 
the conductance as well as other physical quantities can 
not depend on W. Even in the limit of W — ► oo, while 
keeping the (non-interacting) DOS at the Fermi level, t'o, 
constant, the expressions should be finite, since the en- 
ergies involved in transport processes should occur near 
the Fermi energy Ep. Since = £s_^2L ; this scale is 



also ruled out. Indeed, it is possible to express all the 
physical properties using Eg and Sg. Using similar argu- 
ments we can understand the formula for Ec- Since the 
DOS is a static property of the system, determined by 
interaction energies on scales much larger than the local- 
ization length £, it cannot depend on £. The only energy 
scale that we can construct from Eg and Sg that does not 
depend on £ is Eg/Sg = E C - 

Optimizing Eq. with the DOS described above 
yields Mott's law for temperatures above a certain tem- 
perature and ES's law [l| for lower temperatures. The 
crossover temperature in two dimensions is given by : 



TEMPERATURE DEPENDENCE OF THE 
COULOMB GAP 

Determining the value of the finite DOS at the Fermi 
energy due to the effect of temperature is a problem 
which has not been settled yet: numerical investigations 
by [HI seem to be in disagreement with other theoret- 
ical predictions j26[. We hereby present the results of 
the mean-field theory. We calculate the density of states 
(DOS) at the Fermi-energy in the presence Coulomb in- 
teractions and finite temperature. For an infinite system 
at zero-temperature, the Efros-Shklovskii Coulomb gap 
will make the DOS vanish at the Fermi-energy. We know 
that finite-size effects will give rise to a finite DOS, pro- 
portional to ^7=[l2j]. Let us assume we have a system 
large enough so that these finite-size effects will be neg- 
ligible compared to the effects of the finite temperature. 
The finite-temperature mean-field equations are fl3l ]: 



EL 




(4) 



For half filling, Ep = 0, and the positive background 
charge Qb = 0.5. 

We have solved these equations numerically in two di- 
mensions, and found a linear dependence of the DOS at 
Ef as function of temperature. This is consistent with 
ref. (2(|. The results are demonstrated in Fig. [2j Thus 
v{Ef) = aT, and since the Coulomb gap is washed out at 
a temperature of order Ec, we know that a ~ -g 2 -. For 
the data corresponding to Fig. [2] we find that a rts 0.15. 
In this case the ratio m 0.1, which is consistent. We 
can now understand why this temperature dependence 
does not affect the crossover discussed in the previous sec- 
tion: while the finite-temperature contribution is linear 
in T, the optimal hop energy scales as y/T at low temper- 
atures (before the crossover to Mott's regime). One can 
check that even at the crossover temperature the finite- 
temperature correction to the DOS is not important. 



T x = E 3 g/S 2 g. 



(3) 



Note that this is consistent with the physical picture 
given previously, where it was shown that all physical 
properties must depend on Eg and 8g. 

So far we have reviewed simple heuristic arguments and 
discussed Mott's and ES's laws in a non-rigorous way. 
This enabled us to introduce the energy and length scale 
in the problem we study. In the rest of the paper we will 
perform analytical and numerical analysis of the Hartree 
approach, and show that it reproduces the non-rigorous 
results although it includes temperature and correlation 
effects that were omitted in the above heuristic approach. 



MEAN-FIELD CALCULATION OF THE 
CONDUCTANCE 

We calculated the conductance numerically using the 
following scheme: 

1. Find the equilibrium occupation number and energy 
associated with each electronic site within the mean-field 
picture, using Eq. {2J. This takes the interactions into 
account, and the density-of-states obtained will manifest 
the Coulomb gap 

2. Construct the Miller- Abrahams resistor network 



27L 1281 ] . For a non- interacting problem, this is an exact 



method to find the resistance of the system by finding 
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([3]), approximately e - ^- 1 behavior with £ ~ 0.34±0.01 
was observed for the resistivity, close to the Mott value 
in two dimensions of |. Below this temperature C ~ 
0.49 ± 0.02 was observed, in accordance with the Efros- 
Shklvoskii result [1^ . In the next section we give a scal- 
ing argument, which shows that: 
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FIG. 2: Density of states as a function of energy and tem- 
perature. Each graph is an average over 2000 iterations. For 

each, N — 1000, and the ratio — , where T is the temper- 
ed 2 

ature, and n the density, was varied between 0.01 and 0.58. 
The on-site disorder was taken uniform in the interval [—5, 5]. 



that of a certain resistor network. Here, the resistor be- 



tween the sites i and j is given by Ri 



kT 



where 



7? are the equilibrium transition rates , which are calcu- 
lated based on the energies Ei and occupation numbers 
of the sites, see Eq. ([T]). Two sites which are close to each 
other in space and close in energy to the Fermi energy, 
will have a lower value of resistance between them: if 
the sites are far away in space, the overlap integrals will 
vanish exponentially. If the sites are far from the Fermi 
energy, they will be permanently full or empty, and will 
not contribute to the conductance. 

3. Find the resistance of the system. This is done 
by solving for the energies (voltages) at the sites. In 
steady-state, the sum of currents into and out of each 
site must vanish. Expressing these currents in terms of 
the resistor network values and the energies yields a set 
of linear equations, solved numerically. 

Executing the numerical procedure, we found that 
above a certain crossover temperature T x , given by Eq. 



(5) 



This allows us to collapse numerical data obtained for 
different sets of parameters. The crossover between the 
Mott and Efros-Shklovskii regimes is manifested in the 
asymptotics of f(x). Fig. [3] demonstrates the crossover 
between the two regimes, and the data collapse, validat- 
ing the scaling relation. We should emphasize, however, 
that the scaling relation is derived under certain assump- 
tions: we assume that the temperature is low enough 
such that we do not have the trivial scenario of nearest- 
neighbor hopping. At T > W, we shall have nearest- 
neighbor hopping, and Mott's law will not be valid: this 
gives a 'smoothing off' of the Mott exponent, to a tem- 
perature independent regime. Thus, we do not expect 
complete data collapse for the regime where the tempera- 
ture is larger than the disorder. Additionally, at very low 
temperatures, finite-size effects give rise to a finite DOS, 
and thus give rise to deviations from the Efros-Shklovskii 
hopping mechanism. Nevertheless, Fig. [3] shows clearly 
that data collapse is obtained. The least-squares fits in 
the asymptotic (linear) regimes gave the exponents dis- 
cussed previously. 



FORM OF THE CROSSOVER 

The crossover between the Mott and ES regimes was 
discussed in references [l7, 18, 19|, 20j . 

We shall now calculate the form of the crossover within 
our theoretical framework. We shall show that scaling is 
obeyed, for any dimensions, regardless of the details of 
the Coulomb gap. This generalizes previous work, for any 
DOS which depends only on the parameters e 2 and vq. 
As previously explained, this is the generic form for the 
DOS. Notice that in principle the DOS is weakly depen- 
dent on temperature, and therefore this scaling relation 
is only an approximation. Let us proceed to the proof, 
for any dimension d. The basic equation connecting the 
energy E and distance r associated with a hop is obtained 
by requiring that the average number of sites within an 



area of order r 2 and energy lower than E should be of 
order unity. Denoting the DOS as v{E), this gives: 



v(E')dE' ~ 1. (6) 



The second equation comes from optimizing r/ = ^ + 
The optimal r] will determine the conductance of the 
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system via log a ~ —77. Differentiating the first equation 
with respect to r gives 



. , —d dr 



r d+i dE 

Demanding the derivative of 77 with respect to r to vanish 
gives 

f 1 dE 
f + T~dr 

Combining the last two equations yields 



0. 



v(E) 



kd 



[r (E)] d+1 T 



therefore: 
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FIG. 3: The crossover between ES and Mott VRH as the 
temperature is varied. Using the scaling relation of Eq. 
we have collapsed data from different parameter sets onto the 
same master curve. We use N between 500 and 3000, and 
each point represents the average over 10 runs. We have ver- 
ified that logarithmic averaging yields essentially the same 
results. The sites were chosen with uniform probability in a 
s/NXyN square with density n = 1, and the energy scale 

1 — i 

was chosen such that e 77,2 =1. na£ was varied between 0.1 
and 0.2, where £ is the localization length, and the on-site dis- 
order was distributed uniformly in the interval [— |, |]. The 
temperature axis was scaled according to the crossover tem- 
perature of Eq. ([3J. The conductance was scaled according 
to the scaling relation of Eq. (|14|l The lower temperature 
fit has an exponent of 0.49 ± 0.02, consistent with the Efros- 
Shklovskii value of 1/2 while the higher temperature fit has 
an exponent of 0.34 ± 0.01, consistent with the Mott value of 
1/3. The lines are the best fits for the relevant data sets. In 
the Mott regime where the linear fit was taken the conduc- 
tance changes by more than 3 orders of magnitude. In the 
ES regime, the conductance changes by more than 6 orders 
of magnitude. 



To proceed, it is useful to define the function ip(E) 
— — V ^ E " >d d E . We have by assumption ip - 

(ilM)3+T 

and we obtain the equation: 



( j 



[11) 



(8) Since the dimensions of the LHS are [energy /length] d + 1 , 

ip must have the form: ip = g{m)[EivQ e 2 ]' J + T , where 
777. = — — 1 is a dimensionless combination of the pa- 



rameters of ip, and g{m) is a certain function. Then the 
equation takes the form: 



E 



= fJ 



T 



,ET-v$e 2 i, 



(12) 



We have an equation of the form AE = f(BE), with 
A = nj 1 h and B = e d vn(i)i. Both AE and BE 



are dimensionless. The general solution will give us E = 
\4>{^)i for a certain function (f>. Plugging in the values 
of A and B, we obtain 



E 
T 



1 



with T x as defined by Eq. |3|) 



3=T £ 1 



(13) 



and 9 is a non-universal 
function, depending on the form of the DOS. 

A similar dimensional analysis shows that | follows 
an identical scaling law (but with a different function). 
Using the definition of 77, we reach the conclusion that 
the conductance follows the proposed scaling, which can 
be written in terms of the energy scales i5j and E^ as: 



log cr 



f 



d-i 



IT 



(14) 



Notice that for d = 2 we retrieve the crossover temper- 
ature of Eq. ([3]). The analysis did not assume anything 
on the shape of the DOS, other than the assumption that 
it is a function of vq and e 2 : it does not have to contain 
a plateau or a power-law Coulomb gap. Nevertheless, if 
it does contain these features, / must have asymptotics 
which correspond to the Mott and ES VRH: at x 3> 1 
we have / oc i 3 + I i while for 1 « 1 we have / oc I'. 
An advantage of this calculation compared to other the- 
ories is that it allows one to complete the calculation for 
any DOS. Hence, it can be used by experimentalists to 
'reverse-engineer' and find the function / that described 
the DOS from the conductivity measurements. 



SUMMARY 

We have given a consistent Hartree-theory for the con- 
ductance of the electron glass. The long-range nature of 
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the interactions should justify this approximation. In 
a previous work we showed the th eory gives a linear 
Coulomb gap at zero temperature [13| |. in accordance 
with other theories. Since the interactions are treated 
on a mean-field level, we can still use the single-particle 
DOS to characterize the conductance. This explains why 
the Coulomb gap indeed affects the conductance. In- 
deed, by numerically solving the Miller- Abrahams resis- 
tor network for the system, we found a crossover between 
the ES and Mott VRH regimes. Scaling is generically 
obeyed at the crossover, which we have derived analyt- 
ically (neglecting interaction induced correlations), and 
compared to the numerics (which takes them into ac- 
count). In future work we intend to incorporate many- 
electron tunneling into the theory, as well as study the 
out-of-cquilibrium properties of the system. 
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